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Abstract 



We study the formation of vacuum condensates in 2 + 1 dimensional QED in the 
presence of inhomogeneous background magnetic fields. For a large class of magnetic 
fields, the condensate is shown to be proportional to the inhomogeneous magnetic 
field, in the large flux limit. This may be viewed as a local form of the integrated 
degeneracy-fiux relation of Aharonov and Casher. 

1 Introduction 

Parity- and flavor-symmetry breaking aspects of 2 + 1 dimensional QED have been the 
subject of much research in recent years. This subject has applications in planar electron 
systems and also provides a deep analogue of certain features of symmetry breaking in 
3 + 1 dimensional theories relevant for particle physics |1|, 0, Q . An important focus of these 
studies is the question of induced charges and spins, which are also related to induced vacuum 
condensates. Gusynin et al ^ have shown recently that a uniform background magnetic field 
acts as a catalyst for dynamical flavor symmetry breaking in 2 + 1 dimensions. A key part 
of this argument is the appearance of a nonzero vacuum flavor condensate, in the limit of 
zero fermion mass, in the presence of a uniform background magnetic field of strength B: 

< 0|^V|0 > \m-^o = -sign{m)— (1) 

Ztx 

While much can be learned from this constant B field case, in order to include dynamical 
gauge fields it is desirable to have a more complete understanding of this phenomenon for 
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more general background electromagnetic fields. As refers to a uniform B field, it is of 
course consistent with the integrated relation 

/ (fx < 0\ip{x)ilj{x)\0 > \m^o = —sign{m)— [ d?xB{x) = —sign{m)^ (2) 
J zn J 

(where $ is the net magnetic flux), which is essentially Landau's degeneracy- flux relation 
and which was extended to inhomogeneous magnetic fields by Aharonov and Casher |^. 
Much important work has been done exploring the detailed global aspects of this integrated 
result (0), and relating it to mathematical index theorems [|, [TT|, [T^. 

The emphasis of this paper is rather different - here we investigate the extent to which 
(1^) may be viewed as a local relation 

7 1 

< 0\ij{x)ij{x)\0 > \m^o = -sign{m)—B{x) (3) 

/vr 

when the background magnetic field is inhomogeneous. This, and the closely related issue 
of induced spin, have been addressed for the special case of an Aharonov-Bohm flux string 
magnetic field [O. In this paper, we consider the formation of a vacuum condensate in 
the presence of a more general spatially inhomogeneous static background magnetic field. 
We present some illustrative examples in which the condensate may be evaluated explicitly, 
and then we show that for a large class of inhomogeneous magnetic fields the condensate is 
proportional to the magnetic field [just as in @], but only in the large flux limit. 

In Section II we give a brief review of vacuum condensates in 2 + 1 dimensional QED. 
Section III contains two explicit examples of particular inhomogeneous magnetic fields and 
Section IV contains the general asymptotic analysis for radial magnetic fields. Finally, we 
conclude with some brief comments. 

2 Vacuum Condensates in Planar QED 

Consider a parity invariant model of 2 + 1 dimensional quantum electrodynamics with 
fermionic Lagrange density 

Cf = ^ iiT^D^ -m)^ (4) 

Here is a. four- component spinor and the gamma matrices belong to a 4 x 4 reducible 
representation 

where the 2x2 irreducible gamma matrices 7^ are given by 
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These gamma matrices are normalized as {r'^,r'^} = —2g^^l^ where the flat Minkowski 
metric is g^^ = diag{—l, 1, 1). The covariant derivative operator is = df^ — iA^, where 
for notational convenience we have absorbed a factor of "e" into the gauge field A^. This 
model is invariant under the generalized parity transformation [|1^, ^ 
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(7) 

and in the massless limit, m — > 0, has a global U{2) flavor symmetry corresponding to the 
interchange of the two 2x2 irreducible representations. 

We consider static background gauge fields and work in the Weyl {Aq = 0) gauge. Then 
the Dirac equation, [iV^D^ — m)ip = 0, block diagonalizes as 

/ E~m -{Di-iD^) \ 
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E — m 



iIj = 



{Di + iDi) E + m 
E + m 

{Di + tD2) 

which illustrates the fact that this reducible representation model is equivalent to a theory 
describing two species of two-component spinors, one with mass +m and the other with 
mass —m. Without loss of generality, we choose m to be positive, and we also choose the 
net magnetic flux to be positive. 

The upper 2x2 sub-block of the Dirac equation (Q), corresponding to the positive mass 
species, is solved by a two-component spinor 
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when E ^ —m, and where f{x,y) is a solution of the two-dimensional partial differential 
equation 

-{D,-iD2){D,+iD2)f = a'f (10) 

with = E'^ — m^. Thus, there are solutions of positive and negative energy, E = 
=bv^aM-m?. When \E\ m, we can write these solutions (including the appropriate nor- 
malization factors) as 
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where / satisfies the Schrodinger-hke equation (|^). Here, the subscript {1} refers to species 
{1} which corresponds to mass +m, while species {2} has mass —m. The superscripts (±) 
refer to the positive and negative energy solutions. 

The threshold states, with \E\ = m, are special and must be specified separately. Indeed, 
already from (11) we see that for species {1} we can have a positive energy solution with 
but the 1 / ^ 
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m factor excludes a negative energy threshold state of this 
form. By contrast, for species {2} we can have a negative energy solution of this form with 
while the 1 / ^ 



E 



-m, 
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m factor now excludes a positive energy threshold state. 
This imbalance leads to an asymmetry in the spectrum of states, and this asymmetry is the 
ultimate source of the interesting symmetry breaking effects in planar QED. 
The explicit threshold states are 



(0+) 

{1} 



-imt 






V J 



(0-) 
{2} 



+imt 



( \ 



V y 



(12) 



where f^'^\x,y) satisfies the first- order threshold equation^ 

(Di+2D2)/(°^=0 



(13) 



Now expand the fermion field operator in terms of creation and annihilation operators for 
an orthonormal set of positive and negative energy modes from (0) and (|12D- 



^ Jn Jp 



(14) 



^Note that we have excluded the potential threshold states of the form T/jj" j 
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where g^'^^ satisfies (Di — i-D2).9'"-' = 0, because if the solutions to (|l^ ) are normalizable then 



these (7^°-' solutions are not, and vice versa M, I 
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where bn,p and dn,p are fermionic annihilation operators. The label n refers to the eigenvalue 
of the equation (|l^) and hence specifies the energy, while the label p distinguishes between 
degenerate states. Note that both n and p may take discrete and/or continuous values, 
depending on the equations ([10|) and (|13]) respectively. Also, note that the sum over the 



species in (|I^ is understood to include the positive energy threshold states for species {1} 
and the negative energy threshold states for species {2}. 

The vacuum expectation value < 0|^^\l/|0 >=< 0|\l/^r°\I^|0 > is then given by 

< 0|M|0 >= l/fP - 2m ^^^^^ (15) 

where the first term on the RHS only involves threshold states, while the second term involves 
all states with \E\ > m. In the massless limit the second term vanishes and the condensate 
is simply 

< 0\^{x)^{x)\0 > = - i l/f (5^)P (16) 



The minus sign on the RHS is due to the fact that the condensate is a vacuum expectation 
value and so it is a sum over occupied negative energy states; and with m positive, the 
only negative energy threshold states correspond to species {2}, for which the r° sub-block 
is —7°. Changing the sign of m corresponds to interchanging the two species, so one has 
instead +7°. Thus, the condensate should be more precisely written as 

< 0|^'(f)^(f)|0 > U^o = -stgn{m)^\ff\x)\^ (17) 



Note that the condensate is determined entirely by the threshold states, which solve the 
first-order equation (0). We now consider several examples in which this condensate may 
be computed explicitly. 

We begin with the familiar case of a uniform background magnetic field [^]. There is 
still gauge freedom of how we choose to represent the corresponding vector potential. This 
choice of gauge will determine the precise form of the threshold condition (|T3p (as well as the 
eigenvalue equation (|l^) which determines the complete spectrum). In the 'linear gauge', 
with A = {0,Bx), the threshold state equation (|TBp has normalized solutions 

f(^){x,y) = (^^y\^pye-(p~^^)y(^^) (18) 

where the degeneracy label p takes continuous values corresponding to a plane wave in the 
y direction. Thus, it is trivial to evaluate the condensate to be 

< 0|'I'(f)^(f)|0 > U^o = -sign{m)\ - f ^e-^f"^")'/^ = -sign{m)— (19) 

V 7r J 27r 27r 
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In the 'radial gauge', with A = ■j{—y,x), the threshold state equation ( ]T3| ) has normalized 
solutions 



\ 



1 /B\P+^ 
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P^-B\z\y4 



(20) 



where we have defined the complex variable z = x + iy, and the degeneracy label p now takes 
integer values p = 0,1,2,.... Once again, it is trivial to evaluate the condensate to be 



< 0|^(f)^(f)|0 > 



m~>0 



-sign[m)—e 



p=0 



p\ 



-signim)— (21) 



The answer is, of course, the same in each gauge. Also note that in each of these cases, (p!9|) 
and (plD , the condensate is independent of X, cLS IS expected for a uniform B field. We now 
turn to some less trivial cases in which the magnetic field is spatially inhomogeneous. 



3 Inhomogeneous Magnetic Fields: Two Examples 

In this Section we consider two illustrative examples of specific inhomogeneous background 
magnetic fields. The first example is in the 'radial gauge', for which we choose the gauge 
field to be 

A={-dyCt),dA) (22) 

where = (j){r) is some function only of the radial coordinate r. Then the magnetic field is 
radial, B{r) = "^(^^^0); and the (un-normalized but mutually orthogonal) solutions to the 
threshold condition (pTSf ) are 

/f = z^e-'^ (23) 

where p is a non-negative integer. 

We now choose a particular functional form for (p which will permit the explicit normal- 
ization of these states: 

<Pir) = ^log (l + (24) 
The corresponding radial magnetic field is 

^(^) = , ^.x2 (25) 
(1 + *) 

which has a finite net fiux 

$ = i- I d'xB{r) = ^ (26) 
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The constant B represents the maximum value of the magnetic field, and i? is a characteristic 
length scale associated with the spatial variation of the magnetic field. In the infinite fiux 
limit, with i? — oo, this example reduces to the constant B field example. 

With this radial choice ( p4D for 0, the threshold states in (p3D may be normalized, yielding 

1 / 1 \ (P+l)/2 1 ~p 

ei-^y) = ^ 4,) .\ .... (27) 



; 

where f3{u,v) = r{u)r{v) /r{u + v) is the beta function. Having finite fiux, this system 
displays the novel feature that only a finite number of these states are localized and normal- 
izable. Indeed, these states are only localized for p < [$], where [$] is the largest integer 
less than $. The contribution to the condensate from these 'bound' states is 

1 1 (-Y' 

< 0|^(f)^(f)|0 > L^o =-sign(m) ^ V ttt — r (28) 

' ^ ^ ^ ^' ' ^ ^i?2^i + ^)%t^o/5(p+l,$-p-l) ^ ^ 

When $ is an integer we can in fact evaluate this sum exactly, yielding 

< Olt (J)4'(f)|0 > |„„„ = -signim) (l - -L_?_^ (29) 



Note that the condensate has the same form as the magnetic field, but with an overall 
multiplicative factor depending on the net magnetic fiux $. For large $ this factor approaches 
unity, so that 

< 0|^'(f)^(f)|0 > U^o ~ -sign{m)^ (30) 

For our second example, consider a magnetic field that is uniform in one direction (say the 
y-direction), and spatially varying in the x direction. To achieve this type of configuration, 
we choose a convenient 'linear gauge' with A = (0,a(x)). The corresponding magnetic field, 
B{x) = a'{x), is just a function of x, and the (un-normalized but mutually orthogonal) 
solutions to the threshold equation (p!3D are 

/W = eip3/e-rw^)-p) (31) 
The specific choice a{x) = XBtanh{x/X) yields a magnetic field profile 



^The integrated condensate is proportional to the net magnetic flux with both localized and continuum 
states contributing [|[ ^, |ll|, |l2j. Here, for a local analysis of the condensate density we only consider 
the localized bound states, in part because the continuum states contribute at infinity, and also because in 
the large flux limit the magnitude of their contribution is negligible. 
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The corresponding flux is 

$ = — / d^xB = (33) 

271 J TT 

where we have compactified the y-direction with a length L. For this type of magnetic field 
it is, in fact, possible to solve equation (p!OD for the entire spectrum, permitting for example 
the exact evaluation of the effective energy |jl5|. Here, however, for the evaluation of the 
vacuum condensate, we only need the normalized threshold states 



1 



= A -=- , e'Py j-j (34) 

^ ^ LX/p(^B\^+p\,B\^-p\) {2cosh{x/X)f' 



With the ?/-direction compactified, the degeneracy label p takes discrete values p = 
where k is an integer such that 

W < 4 (35) 



in order for these states to decay at infinity. We can therefore perform the sum in (p!7D, 
yielding 



\ BLXl 



iT.^xw^M , 2sign(m) ink{x/L) 

<0^f^f 0>L^0= on.2 y ^ X ^ 



LA {2cosh{x/\)f''^' (3{B\^ + 27rfc^, EA^ - 27rA;^) 

(36) 

As L — >■ oo we can replace the sum over k by an integral and find 

(37) 

It is straightforward to plot this condensate for various values of the dimensionless com- 
bination BX^. One finds that the condensate has the same general 'bell-shaped' form as 
—sign{m)B{x) /27r, but that the correspondence is not exact. Nevertheless, for large BX^ 
(which corresponds to large flux) we can use Stirling's formula to make an asymptotic ex- 
pansion of the inverse beta function in the integrand of (|37D to obtain 



B \ BX^ 1 

< 0|^(f)^'(f)|0 > U^o ~ -sign{m) — \j—-- _ ,^.^2BX^ (3^ 



27r 



TT {cosh{x/X)y 



j dtVl — t'^exp 



BX^{2t{^)-log{l-t^)-tlog^^^^ 



.XJ ' "Vl-t 

The remaining integral over t is suited for an asymptotic expansion, for large BX^, using 
Laplace's method [0. For an integral of the form 

I{N) = f ds'^{s)exp{Nn{s)) (39) 
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the large leading asymptotic behavior is dominated by a critical value Sc at which the 
exponent function Q[s) has a maximum, and is given by 



/ 271 

^^^) ~ ^ -NQ"{s,) "^^'^^ ''''^ (N^i^c)) (40) 

Applying Laplace's method to the t integral in (|39|) , for which the critical point is at tc = 
tanh{x/X), we find the condensate to be asymptotically proportional to the inhomogeneous 
magnetic field: 

B 1 B(x) 

< 0|*(f)^'(f)|0 > l^^o ~ -stgn{m)-— . = -stgn{m)^^ (41) 

^TT [cosh{x / X)) 2tc 

4 Asymptotic Analysis for General B{r) 

At first sight, one might think that the asymptotic proportionality between the condenstate 
and the inhomogeneous magnetic field found in (|30| ) and ( |^ ) is due to the special form of 
the inhomogeneous magnetic field chosen in these examples. For example, in each of these 
cases the normalization factors may be computed exactly and are given by beta functions. 
In general it is not possible to compute the normalization factors in closed form. However, 
we show in this Section that this asymptotic analysis applies to very general inhomogeneous 
magnetic fields B{r) and B{x). As the analysis is very similar in the two cases, we concentrate 
on the radial case. In the conclusion we make some comments concerning the general B{x, y) 
case. 

We choose the gauge field in the radial gauge (^) with = 0(r). It is convenient to 
write 

BR^ /r^\ 

= -^h (^^ j (42) 

where R is some characteristic length scale, and B together with the dimensionless function 
h are chosen so that the overall normalization gives net flux $ = "^^^ special cases 

^(0 = i ^-iid h{^) = log{l + ^) have been considered in the preceeding Sections. Then, as 
before, the threshold states are given by (PBf). Note that these are automatically mutually 
orthogonal for any (j){r) by virtue of the angular integration. The vacuum condensate is 

('^//?)2Pg-<I'h{r2/R2) 

< 0|^'(x)^(x)|0 > U^o = -sign{m) J2 1^2 (^3) 

p=0 ^> 

where the normalization factors A^^ are given by 



p 

v#J \r) 



N' = nR' I 'd U-] e-'^HrVR^) (44) 



9 



For large flux, the sum over p in (^3]) may be replaced by an integral from to $, which 
may be re-expressed in terms of the rescaled variable t = p/$ as 



< 0|^(f)^'(f)|0 > U^o = -stgn{m)(!> dt i j (45) 

with normalization factors _ 

poo 

N^ = nR^ rf^e*(*'°3«-'^«« (46) 
Jo 



V.2 



where ^ = 

To evaluate the asymptotic form of the condensate (^) for large flux $ we first need the 
asymptotic form of the normalization factors for large $. This can also be done by Laplace's 
method, with the dominant contribution coming from the maximum of the exponent function 

X{0=tlog^-m (47) 

This defines a critical point = C,c{t), as a function of the parameter t, via the implicit 
relation 

^h'iO = t (48) 

Using (|40|) we see that 

Inserting this into the integrand in (|45|) we obtain 



< 0|f^(x)^(f)|0 > [m^o ~ -sign{m) — \ — + /i"(^^(t))e*[«°9«/«=W)-^(«)+'^(«=W)l 



27r V 27r Jo V ic{t) 

(50) 

The remaining t integral may also be expanded asymptotically using Laplace's method, with 
exponent function 

fi(t) = tlogi - tlogUt) - HO + HUt)) (51) 

This leads to a maximum at the critical point tc = tc(0 a function of ^ = r"^ / E?, defined 
by the relation ^ = C,c{tc)- Applying the inverse function theorem (see comments below 
equation (|51D), we find 

tc = ^h'iO (52) 
and ^ 

~ ^"^''^ = mo+mo) ^^^^ 



10 



Combining all these pieces, several remarkable cancellations occur, and the leading asymp- 
totic form of (^) is simply given by 

< 0\^{x)^{x)\0 > U^o sign{m)^ {h'{0+^h"{0) 

ZTC 

Bir) 

= —sign{m) (54) 

27r 

This result depends on subtle cancellations between the asymptotic expansion of the (in- 
verse of the) normalization integral (which is an integral over r or equivalently ^) and the 
asymptotic expansion of the sum over the threshold states (which becomes an integral over 
t). These cancellations rely on the use of the inverse function theorem which assumes that 
the function Ch'{^) is one-to-one (see (|48|)). This means that its derivative, (^/i'(0)') has a 
fixed sign, which moreover must be positive in order for the critical point to be a maximum. 
Since (^/i'(0)' is just the magnetic field, we see that our result applies over a region of space 
in which the inhomogeneous magnetic field is positive. This is the relevant physical set-up 
and is consistent with our intuitive expectation - in a region with positive magnetic field, the 
spin density tends to align with the magnetic field; moreover, if the flux in this region is very 
large, the local spin density will actually be approximately proportional to the local inhomo- 
geneous magnetic field. The proportionality factor is such that when this result is integrated 
over the region, we regain the Aharonov-Casher relation |^ between the integrated spin 
and the net magnetic flux. 



5 Concluding Remarks 

In this paper we have studied the vacuum condensate of parity invariant 2 + 1 dimensional 
QED in the presence of inhomogeneous magnetic fields. After presenting two explicit illus- 
trative examples in which the condensate may be evaluated in detail, we showed that in the 
limit of large flux the condensate is locally proportional to the inhomogeneous magnetic field: 

{—sign{m)^^ 
(55) 
/ ^B(x) 
-sign[m)^ 

for general physically relevant magnetic fields which either depend only on the radial coordi- 
nate or only on one of the Cartesian coordinates. These relations represent a local analogue of 
the integrated Aharonov-Casher relation. To generalize this result to general static magnetic 
fields B{x,y), we note that the background vector potential may be represented as 

A = {-dy(f){x, y), d^(f){x, y)) (56) 
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where 



(57) 



Then the threshold solutions are simply 



(58) 



where Fp{z) is some holomorphic function. However, to proceed one must construct an 
orthogonal basis of such holomorphic functions. When (p = (f){r) this may be achieved by the 
choice Fp{z) = z^, as in ([23|), and when = 0(x) by the choice Fp{z) = e^^, as in (pT)). For a 
general (p = (f){x, y) the natural choice of orthogonal basis is not so clear (a Gramm-Schmidt 
orthogonalization is too clumsy for the subsequent summation). Further clarification of this 
issue should lead to interesting insights into the properties of vacuum condensates of 2 + 1 
dimensional QED with djTiamical gauge fields. 
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